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Floating car data of car-following behavior in cities were 
compared to existing microsimulation models, after their pa- 
rameters had been calibrated to the experimental data. With 
these parameter values, additional simulations have been car- 
ried out, e.g. of a moving car which approaches a stopped 
car. It turned out that, in order to manage such kinds of sit- 
uations without producing accidents, improved traffic models 
are needed. Good results have been obtained with the pro- 
posed generalized force model. 
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I. INTRODUCTION 

During the last five years, gas- kinetic fluid- 
dynamic and other models have been developed, 
aiming at an understanding of stop-and-go traffic. The 
topic is related to the fields of non- linear dynamics 
phase transitions , and stochastic processes |^ . In ad- 
dition, microscopic traffic models were proposed for the 
description of interacting driver- vehicle units. They can 
be classified into cellular automata models Q , which are 
discrete in space and time, and continuous models . 
The latter are required for detail studies of car-following 
behavior and traffic instabilities, which are necessary for 
an investigation of the consequences of technical opti- 
mization measures (e.g. of autopilots for an automatic 
control of vehicle acceleration and braking). 

Therefore, a research group of the Bosch GmbH has 
recently recorded follow-the-leader data by means of a 
floating car a which measured the vehicle speed Va, the 
netto distance Sa to the car in front, the acceleration aa 
of it, and the relative velocity Ava- By a correlation anal- 
ysis it was demonstrated that, among all possible combi- 
nations of subsets of these four quantities, Sq., Ava, and 
Va are the most significant variables for the description 
of vehicle dynamics Q. In Sec. Ill, we will find plausible 
reasons for this. 



The follow-the-leader data, if plotted in the Sa-Ava- 
plane, show the characteristic oscillation of vehicle mo- 
tion around states with relative velocity zero (cf. Fig. ||), 
which was already reported by Hoefs ||^ . Except for the 
previously mentioned significance analysis, the data were 
also used for calibrating existing microsimulation models. 
With the resulting optimal sets of parameter values, the 
models were simulated for the observed situation. That 
is, the first vehicle was moved according to its measured 
velocity, and the following vehicle was simulated accord- 



ing to the respective model under consideration, starting 
with the the same initial velocity and distance like the 
floating car. The average relative quadratic deviation D 
between the simulated and actually measured distance 
(cf. Eq. (p^) was used as a measure for the goodness of 
fit of the respective model Q . 

In order to get improved results, we have developed 
a generalized force model, in which each term and each 
parameter has a clear meaning. Moreover, by parameter 
calibration it turns out that all model parameters have 
the right order of magnitude. Therefore, it can be eas- 
ily said, how the parameter values will look like, if the 
speed limit, the acceleration capability, the average vehi- 
cle length, the visibility, or the reaction time is modified 
(e.g. due to technical measures). In addition, this model 
reaches a better goodness of fit at a reduced number of 
model parameters than the previous models. Finally, 
the generalized force model manages to cope success- 
fully with particular situations like vehicles approaching 
standing cars, in which other models produce accidents. 



II. DISCUSSION OF PREVIOUS MODELS 

The first microscopic traffic models were developed in 
the 1960ies. Many of them are special cases of the follow- 
the-leader model proposed by Gazis, Herman, and Roth- 
ery This assumes that the dynamics of a vehicle a 
with velocity Va (t) at place Xa (t) is given by the equation 
of motion 



dxa (t) 
dt 



Va{t) 



and the acceleration equation 
dva{t + T) 



dt 



K{t + T)[VaMt)-Va{t)]. 



(1) 



(2) 



According to this, a driver adapts to the velocity Va-i(t) 
of the car in front, but this is delayed by the adaptation 
time T « 1.3 s. The deceleration is proportional to the 
relative velocity 



AVa = - Va-1 



(3) 



where the proportionality factor k reflects the sensitivity 
to the stimulus Ava- The sensitivity was assumed to 
depend on the vehicle velocity and on the brutto distance 



in the following way: 



(4) 



1 



K{t + T) ~ Kq 



[Va{t + T)Y 



(5) 



This choice allowed to fit all equilibrium velocity-density 
relations of the form 



Vc{p) = Vq 



1 - 



Pn 



1/(1-™) 



(6) 



by appropriate specification of the exponents I and m (vq 
= maximum velocity, p = spatial vehicle density, Pmax 
= maximum vehicle density). The best fit was reached 
for fractional exponents m w 0.8 and I « 2.8, so that 
this model has no obvious interpretation. Apart from 
that, the model does not allow to distinguish drivers with 
different preferred velocity, and it cannot describe the 
acceleration of a single vehicle correctly. 

Only a few years ago, Bando, Sugiyama et al. proposed 
a very charming microscopic traffic model. Despite its 
simplicity and its few parameters, their optimal velocity 
model (OVM) describes many properties of real traffic 
flows and is easily interpretable. It is based on the 
acceleration equation 



dvg jt) 
dt 



= K[V{Sa) - Va(t)] , 



(7) 



SO that the vehicles adapt to a distance-dependent opti- 
mal velocity 

Visa) ^Vi+V2 tanh(Cis„ - C2) (8) 

with a certain relaxation time r — Here, 

Sq — Xq^ — \ ^Q — 1 Sa la — 1 (9) 

denotes the netto distance, where la means the length of 
vehicle a. Like the follow-the-leader models, the optimal 
velocity model is able to describe the formation of stop- 
and-go waves and emergent traffic jams, but it overcomes 
the afore mentioned problems. 

We carried out a calibration of the optimal velocity 
model with respect to the empirical follow-the-leader 
data, which we obtained from Bosch. The optimization 
procedure based on the evolutionary Boltzmann strat- 
egy |lo|] , and the optimization criterion was the average 
relative quadratic deviation 



D 



1 ^ 
N ^ 



t=i 



Sait)~S^{t) 



(10) 



of the simulated distance Sa{t) from the measured ve- 
hicle distance s™(i). The resulting optimal parameter 
values for city traffic in Stuttgart are k = 0.85 s~^, 
Vi = 6.75 m/s, V2 = 7.91 m/s, Ci = 0.13 m"!, and 
C2 = 1.57. 

A comparison with the data shows that the extremely 
short relaxation time r = 1/k = 1.17s results in too high 
values of acceleration, which leads to an overshooting of 



the vehicle velocity (cf. Fig. ^) . The unrealistically high 
accelerations also become obvious in Figure |[ since em- 
pirical accelerations are limited to 4m/s^ (cf. Fig. ||c). A 
similar problem occurs with the deceleration behavior, if 
a standing car (e.g. at the end of a traffic jam or in front 
of a red traffic light) is approached from a large distance 
by an initially freely moving car. It turns out that the 
moving vehicle reacts too late to the vehicle at rest. The 
values of deceleration are unrealistic large, but still not 
sufficient to avoid an accident (cf. Fig. 

These problems are solved by the T3 model 

dVa 1 + biVa + 62SQ -I- b^VaSa + b^V 

dt Co + CiVa + C2Sa + C^VaSa + C^Va-l + C^VaVa^l 

(11) 

proposed by Bosch b^ and Ck are model parameters. 
The regression model (pT|), which is based on a ratio- 
nal function, describes all aspects of vehicle dynamics in 
cities realistically (cf. Figs. ||-^, but at the cost of addi- 
tional parameters. Whereas the optimal velocity model 
needs only 5 model parameters, the T3 model contains 
11 parameters. If the model equations are scaled to di- 
mensionless equations (by scaling space or velocity and 
time by characteristic model quantities), the number of 
parameters is reduced by 2. 

In the following section, we will propose an alternative 
model which reaches about the same goodness of fit as 
the T3 model, but with a considerably smaller number 
of parameters. 



III. THE GENERALIZED FORCE MODEL 

Motivated by the success of so-called social force mod- 
els in the description of behavioral changes ||ll|,|l[, es- 
pecially of pedestrian dynamics |l2|Jll], we developed a 
related model for the dynamics of interacting vehicles. 
In setting up an equation of motion by specifying the ef- 
fective acceleration and deceleration forces, the approach 
is analogous to the molecular-dynamic method which is 
used for the simulation of many-particle systems [ p^ , 
e.g. of driven granular media Jl^ . 

Besides of various methodological similarities in the 
theoretical treatment of traffic and granular flows, there 
are also phenomenological analogies like the formation of 
density waves . In both cases, the interactions are dis- 
sipative, i.e. they do not conserve kinetic energy. How- 
ever, there are also differences. For granular media, the 
interaction forces are short-ranged and belong to collision 
processes, where particles touch and temporarily deform 
each other. Vehicle interactions are long-ranged and cor- 
respond to deceleration maneuvers. They are usually not 
related to collisions (i.e. accidents), since the drivers try 
to keep a safe distance s from each other which can be 
considerably larger than the vehicle length (cf. Eq. (p^). 
Therefore, the effective space requirements of vehicles are 



2 



much larger than their actual size. Moreover, vehicular 
interactions do not conserve momentum (in contrast to 
granular ones). 

Another difference between granular and traffic dy- 
namics is, that the laws of granular interactions are very 
well known, pretty much like the basic laws of physics, 
whereas the laws of driver- vehicle dynamics (if they exist 
at all) are still to be established. The particular chal- 
lenge of modeling vehicle dynamics is its dependence on 
factors like perceptions, psychological motivations and 
reactions, or social behaviors. Thus, in contrast to phys- 
ical processes, driver behavior cannot be expected to be 
describeable by a few natural constants. 

According to the social force concept, the amount and 
direction of a behavioral change (here: the temporal 
change of velocity, i.e. the acceleration) is given by a sum 
of generalized forces. These reflect the different motiva- 
tions which an individual feels at the same time, e.g. in 
response to their respective environment. Since these 
forces do not fulfill Newton's laws like actio = reactio, 
they are called generalized forces. Alternatively, they are 
named behavioral or social forces, since they mostly corre- 
spond to social interactions. The success of this approach 
in describing traffic dynamics is based on the fact that 
driver reactions to typical traffic situations are more or 
less automatic and determined by the optimal behavioral 
strategy (which is the results of an initial learning pro- 
cess). A detained motivation, description, and discussion 
of the social force concept is given in Refs. pT|,p],p^ . 

The driver behavior is mainly given by the motiva- 
tion to reach a certain desired velocity (which will be 
reflected by an acceleration force /") and by the motiva- 
tion to keep a safe distance to other cars (3 (which will 
be described by repulsive interaction forces fa.p)'- 



fa.A^a^Va^Xil^Vp) + S,a{t) ■ (12) 



The fluctuating force S,a{t) may be used to include indi- 
vidual variations of driver behavior, but in our present 
investigations it was set to zero. If we assume that the ac- 
celeration force is proportional to the difference between 
desired and actual velocity and suppose that the most 
important interaction concerns the car (a — 1) in front, 
we end up with 

dVa vl - Vg ^ f , \ f,o\ 

—7—^ h Ja.a-l(Xa,Va;Xa-l,Va-l) ■ (lo) 

dt Tq, 

The acceleration time is a third of the time which 
a freely accelerating vehicle needs to reach 95% of the 
desired velocity. 

Now, we have to specify the interaction force fa.a-i- 
For 



fa,a — l — 



(14) 



and Ta = 1/k, we would obtain the optimal velocity 
model, again. We extend this relation by a complemen- 
tary term which should guarantee early enough and suf- 
flcient braking in cases of large relative velocities Ava- 
This term should increase with growing velocity differ- 
ence Ava, but it should be only effective, if the velocity 
of the following vehicle is larger than that of the leading 
vehicle, i.e. if the Heaviside function Q{Ava) is equal to 
1. Moreover, the additional deceleration term should in- 
crease with decreasing distance s^, but vanish for large 
distance Sq — > oo. The braking time belonging to 
this term should be smaller than t^, since deceleration 
capabilities of vehicles are greater than acceleration ca- 
pabilities. We chose the following formula which meets 
the above conditions: 



/a, a- 



Vis^) 



AVgeiAVg) _[ 



(15) 



This formula takes into account that vehicles prefer to 
keep a certain velocity-dependent safe distance 



s{Va) = da+ TaVa , 



(16) 



where da is the minimal vehicle distance and Ta is the 
safe time headway (i.e. about the reaction time). R'^ can 
be interpreted as range of the braking interaction. 

We can further reduce the number of parameters (and 
the numerical effort), if we replace the previous V{sa)- 
function (||) by 



Va{Sa,Va) = l^" { 1 



-lSc-sM]/R.a 



} 



(17) 



In case of identical model parameters of all vehicles, the 
corresponding equilibrium velocity-distance relation re- 
sults from the implicit condition Va = Va{sa,Va) and is 
depicted in Figure ^. 

The traffic model defined by equations (|l^), ([Tsl), and 
( p7| ) will, in the following, be called the generalized force 
model of traffic dynamics (GFM). Since all its seven pa- 
rameters have a clear and measurable meaning, they 
should have the right order of magnitude. A calibra- 
tion with respect to the follow-the-leader data shows 
that this is indeed the case. We found the following op- 
timal parameter values: w° — 16.98 m/s, Ta = 2.45 s, 
da = 1.38m, Ta = 0.74s, < = 0.77s, Ra = 5.59m, and 
R'a = 98.78 m. Now, the acceleration time Ta is more 
than twice as large as in the optimal velocity model, the 
braking time is smaller than Ta, as demanded, and the 
reaction time Ta is also realistic. Note that the range Ra 
of the acceleration interaction is much shorter than the 
range R'a of the deceleration interaction. This is not only 
sensible, it is also the reason for the astonishingly good 
agreement with the empirical data. Table | compares the 
minimal values of the average relative quadratic deviation 
D that could be reached for the different discussed traf- 
fic models by evolutionary parameter optimization pO[. 
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(The advantage of the appUed Boltzmann strategy is that 
this particular gradient method escapes local minima by 
means of a fluctuation term with eventually decreasing 
variance.) 

It turns out that the optimal velocity model is consid- 
erably improved by the T3 model. This is not surprising, 
since the goodness of fit should increase with the num- 
ber of model parameters. Nevertheless, the generalized 
force model reaches the best agreement with the data, 
although it includes only two third of the number of pa- 
rameters of the T3 model. The simulation results for 
the generalized force model are depicted in Figure |^. Fi- 
nally, the representation of the relative vehicle movement 
in the Aua-Sa-plane shows the expected oscillatory char- 
acter of foUow-the- leader behavior, which can cause the 
development of stop-and-go traffic (cf. Fig. ^). 

In comparison with physical models, the generalized 
force model still seems to contain a lot of parameters. 
However, let us discuss this in more detail for the previ- 
ously mentioned molecular-dynamic models of granular 
media. If we want to describe interactions of smooth, 
inelastic, and spherical grains, only, we need to know 
the particle size and the normal restitution coefRcient 
related to translational energy dissipation. In cases of 
rough spheres, we require two additional parameters, a 
tangential restitution coefficient and a friction coefficient 
[ p6t . Moreover, if the grains are non-spherical, the situa- 
tion becomes even more complicated. In conclusion, the 
small number of parameters occuring in physical models 
are often a result of simplifications and idealizations. 

Finally, let us check the plausibility of the generalized 
force model. In order to do this, we rewrite the model in 
the form 



dVg jt) _ Vg-lit) - Vajt) 

dt r" 



(23) 



with 



and 



where 



dVg _ V*{Sa,Va, AVg) - Va{t) 

dt ~ T* 



1 1 eiAvg) 



< = exp{[s„ - s{Va)]/R'a} ■ 



(18) 



(19) 



(20) 



(21) 



For small velocity differences Ava or large distances Sq, 
we find 



dVajt) ^ Va{Sa,Va) - Vajt) 
dt Try 



(22) 



so that vehicles try to approach the optimal velocity Va 
with a relaxation time Tq. This corresponds to the opti- 
mal velocity model, but with a velocity-dependent func- 
tion Va{sa, Va). If a veliiclc is faster than the leading ve- 
hicle (i.e. Ava > 0) and its distance is sufficiently small, 
we have 



which coincides with a car-following model in which the 
following vehicle adapts to the velocity of the vehicle in 
front. According to the formula for r^', the decelera- 
tion is the stronger, the closer the vehicles come to each 
other. Therefore, the limiting cases of the generalized 
force model behave very reasonably. 



IV. SUMMARY AND DISCUSSION 

We have calibrated several microscopic traffic models 
to city traffic, compared them with empirical foUow-the- 
leader data, and investigated their properties. It turned 
out that one model showed too large accelerations and 
decelerations, but nevertheless caused accidents in cer- 
tain situations. Another model was a regression model, 
so that the meaning of the model and its parameters was 
not clear. Therefore, we developed the generalized force 
model, which reached the best agreement with the em- 
pirical data, although it had only two parameters more 
than the optimal velocity model and four parameters less 
than the T3 model. Another advantage of the general- 
ized force model is that all its parameters are easily in- 
terpretable and have the expected order of magnitude. 
Therefore, it can be immediately said, how the param- 
eters will differ between fast cars, slow cars, and trucks 
(the latter being characterized by small , but large 
and r^). It can also be predicted what happens, if the 
speed limit (i.e. v^) is changed, if the weather conditions 
are bad (greater r^, but smaller and R'^), if roads 
would be used by vehicles with smaller length la, if the 
reaction time Ta could be reduced (by means of technical 
measures like an autopilot), etc. For these reasons, the 
generalized force model is an ideal tool for carrying out 
detail studies of traffic flow, as well as for developing and 
testing traffic optimization measures. 

The simulation of large vehicle numbers is completely 
analogous to molecular-dynamic simulation studies of 
many-particle systems, e.g. of granular flows. The pa- 
rameters of the different driver- vehicle units are specified 
individually (a-dependent), then. In this case, one would 
specify typical parameter values of fast cars, slow cars, 
and trucks, and their respective percentages. Alterna- 
tively, one could introduce a distribution of each param- 
eter (e.g. a Gaussian one) around a typical value. Then, 
the simulation is started with the initial and boundary 
conditions of interest. Of course, the model can be also 
extended to a multi-lane model with lane-changing and 
overtaking maneuvers [l|Jl7||. This is a topic of current 
research. 
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FIG. 1. The foUow-the-leader data show the oscillatory na- 
ture of the relative motion of vehicles. 
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FIG. 2. The time-dependent velocity Vci{t) (a), distance 
Sa{t) (b), and acceleration dva{t)/dt (c) according to the op- 
timal velocity model (OVM), the general force model (GFM), 
and the T3 model in comparison with foUow-the-leader data 
of city traffic. According to (a), most of the models compare 
well to the measured velocities. Only the OVM shows a signif- 
icant overshooting, indicating too large accelerations. Fitting 
of the vehicle distances is a much harder task, as shown in 
(b), but the T3 model and the GFM perform well. In (c), one 
can see that the empirical accelerations and decelerations are 
usually limited to the range between — 3m/s^ and +4m/s^, 
which is met by the GFM. 

Note that the vehicles had to stop three times due to red 
traffic lights (during the periods between 169.9 and 184 sec- 
onds, 233.5 and 253.5 seconds, and after t = 288 s). At time 
t = 143.8 s, the vehicle in front was turning right, so that the 
second car was following another vehicle, afterwards. 
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FIG. 3. Acceleration of an unobstructed vehicle and of 
a following vehicle according to the optimal velocity model 
(OVM) and the T3 model. Initially, both vehicles are at rest. 
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FIG. 4. Time- dependent velocity Va{t) (a) and accelera- 
tion dva (t) /dt (b) for a vehicle which approaches a standing 
vehicle according to the different discussed simulation mod- 
els. The optimal velocity model produces an accident at time 
t = 34.7 s. 
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FIG. 5. Velocity-distance relations of the different traffic 
models in the stationary case, if all vehicles have identical 
parameters. 
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FIG. 6. The simulation of the foUow-the-leader behavior 
according to the generalized force model shows the oscillatory 
nature of the relative motion of vehicles. The above result is 
in good agreement with the empirical findings depicted in 
Fig. I 



TABLE I. Minimal values of the average relative deviation 
D between empirical data and simulation results that were 
reached for the different traffic models by evolutionary pa- 
rameter optimization. 



Model 


OVM 


T3 


GFM 


D 


0.0586 


0.0354 


0.0316 



6 



